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Computation of Thermally Perfect Properties of Oblique Shock Waves 


Summary 

A set of compressible flow relations describ- 
ing flow properties across oblique shock waves, 
derived for a thermally perfect, calorically 
imperfect gas, is applied within the existing 
thermally perfect gas (TPG) computer code. The 
relations are based upon a value of c p expressed 
as a polynomial function of temperature. The 
updated code produces tables of compressible 
flow properties of oblique shock waves, as well as 
the original properties of normal shock waves 
and basic isentropic flow, in a format similar to 
the tables for normal shock waves found in 
NACA Rep. 1135. The code results are validated 
in both the calorically perfect and the calorically 
imperfect, thermally perfect temperature 
regimes through comparisons with the theoreti- 
cal methods of NACA Rep. 1135, and with a 
state-of-the-art computational fluid dynamics 
code. The advantages of the TPG code for 
oblique shock wave calculations, as well as for 
the properties of isentropic flow and normal 
shock waves, are its ease of use, and its applica- 
bility to any type of gas (monatomic, diatomic, 
triatomic, polyatomic, or any specified mixture 
thereof). 

The TPG code computes properties of 
oblique shock waves for given shock wave angles, 
or for given flow deflection (wedge) angles, 
including both strong and weak shock waves, 
including the special case of a normal shock. 
The code also allows calculation of a sequence of 
shocks where each successive shock wave is 
dependent upon the flow conditions from the 
preceding one. Both tabular output and output 
for plotting and post-processing are available. 
The user effort required is minimal, consisting of 
simple answers to 14 interactive questions (or 
less) per table generated. 

Introduction 

The traditional computation of one-dimen- 
sional (1-D) isentropic compressible flow proper- 
ties, and properties across normal and oblique 
shock waves, has been performed with calori- 
cally perfect gas equations such as those found 
in NACA Rep. 1135 (ref. 1). When the gas of 


interest is air and all shock waves present are 
normal to the flow direction, the tables of com- 
pressible flow values provided in NACA Rep. 
1135 are often used. For shock waves in air that 
are oblique to the flow direction, NACA Rep. 
1135 contains charts of certain flow properties. 
These tables and charts were generated from the 
calorically perfect gas equations with a value of 
1.40 for the ratio of specific heats y. The applica- 
tion of these equations, tables, and charts is lim- 
ited to that range of temperatures for which the 
calorically perfect gas assumption is valid. How- 
ever, many aeronautical engineering calcula- 
tions extend beyond the temperature limits of 
the calorically perfect gas assumption, and the 
application of the tables or equations of NACA 
Rep. 1135 can result in significant errors. These 
errors can be greatly reduced by the assumption 
of a thermally perfect, calorically imperfect gas 
in the development of the compressible flow rela- 
tions. (For simplicity within this paper, the term 
thermally perfect will be used to denote a ther- 
mally perfect, calorically imperfect gas.) Previ- 
ous papers (references 2 and 3) described a 
computer code, and the underlying mathemati- 
cal formulation, which implements 1-D isen- 
tropic compressible flow and normal shock wave 
relations derived for a thermally perfect gas. 
The current paper presents an enhanced com- 
puter code, and the corresponding mathematical 
derivation, for the computation of the oblique 
shock wave relations based upon the assumption 
of a thermally perfect gas. 

A calorically perfect gas is by definition a 
gas for which the values of specific heat at con- 
stant pressure c p and specific heat at constant 
volume c v are constants. Reference 7, as well as 
many compressible flow textbooks, derive and 
summarize calorically perfect compressible flow 
relations based upon this definition. The accu- 
racy of these equations is only as good as the 
assumption of a constant c p (and therefore a con- 
stant y). For any non-monatomic gas, the value 
of c p actually varies with temperature and can 
be approximated as a constant for only a rela- 
tively narrow temperature range. As the tem- 
perature increases, the c p value begins to 
increase appreciably due to the excitation of the 
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vibrational energy of the molecules. For air, this 
phenomenon begins around 450 to 500 K. The 
variation of c p with temperature (and only with 
temperature) continues up to temperatures at 
which dissociation begins to occur, approxi- 
mately 1500 K for air. Thus, air is thermally 
perfect over the range of 450 to 1500 K and 
application of the calorically perfect relations 
can result in substantial errors. A similar range 
can be defined for other gases over which the gas 
is calorically imperfect, but still thermally per- 
fect. At yet higher temperatures c p becomes a 
function of both temperature and pressure, and 
the gas is no longer considered thermally per- 
fect. 

NACA Rep. 1135 does present one method 
for computation of the 1-D compressible flow 
properties of a thermally perfect gas. In this 
approach, the variation of heat capacity due to 
the contribution from the vibrational energy 
mode of the molecule is determined from quan- 
tum mechanical considerations by the assump- 
tion of a simple harmonic vibrator model of a 
diatomic molecule (ref. 4). With this assumption 
the vibrational contribution to the heat capacity 
of a diatomic gas takes the form of an exponen- 
tial equation in terms of static temperature and 
a single constant 0. The complete set of ther- 
mally perfect compressible flow relations based 
upon this model is listed in NACA Rep. 1135 
“Imperfect Gas Effects” Tables of these ther- 
mally perfect gas properties are not provided 
because each value of total temperature T t 
would yield a unique table of gas properties. 
Instead, NACA Rep. 1135 provides charts of the 
1-D thermally perfect air isentropic and normal 
shock properties normalized by the calorically 
perfect air values and plotted versus Mach num- 
ber for select values of total temperature. For 
oblique shock waves, an even more limited set of 
charts is presented for shock wave angle a, 
downstream Mach number M 2 , and pressure 
coefficient as functions of deflection angle 8 for 
four static temperatures at two total tempera- 
tures. This approach can only be applied to 
diatomic gases (e.g., N 2 , 0 2 , and H 2 ) because the 
functional form which describes the variation of 
heat capacity with temperature accounts only 
for the harmonic contribution. 


Because the imperfect gas method of NACA 
Rep. 1135 is applicable only to diatomic gases, a 
different method of computing the 1-D isentropic 
and normal shock flow properties of a thermally 
perfect gas was developed and described in refer- 
ences 2 and 3 . The current paper extends that 
methodology to the computation of flow proper- 
ties across oblique shock waves. The method uti- 
lizes a polynomial curve fit of c p versus 
temperature to describe the variation of heat 
capacity for a gas. The data required to generate 
this curve fit for a given gas can be found in tab- 
ulated form in published sources such as the 
NBS “Tables of Thermal Properties of Gases” 
(ref. 5) and the “JANAF Thermochemical Tables” 
(ref. 6). Actual coefficients for specific types of 
polynomial curve fits are also published in NASP 
TM 1107 (ref. 7), NASA SP-3001 (ref. 8 ), and 
NASA TP-3287 (ref. 9). Use of these curve fits 
based upon tables of standard thermodynamic 
properties of gases enables the application of 
this method to any type of gas: monatomic, 
diatomic, and polyatomic (e.g., H 2 0, C0 2 , and 
CF 4 ) gases or mixtures thereof, provided a data 
set of c p versus temperature exists for the indi- 
vidual gas(es) of interest. 

In this report, a set of thermally perfect gas 
equations, derived for the specific heat as a poly- 
nomial function of temperature, is applied to the 
calculation of flow properties across oblique 
shock waves. This set of equations was coded 
into a previously developed computer program 
referred to as the Thermally Perfect Gas (TPG) 
code. The new oblique shock wave capability 
was added as an optional output to the isen- 
tropic flow and normal shock wave tables of the 
original code. All of the output tables of the TPG 
code are structured to resemble the tables of 
compressible flow properties that appear in 
NACA Rep. 1135, but can be computed for arbi- 
trary gases at arbitrary Mach numbers and/or 
static temperatures, for given total tempera- 
tures. All properties are output in tabular form, 
thus eliminating the need for graphical interpo- 
lation from charts. As in the original TPG code, 
the biggest advantage is the validity in the ther- 
mally perfect temperature regime as well as in 
the calorically perfect regime, and its applicabil- 
ity to any type of gas (monatomic, diatomic, tri- 
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atomic, polyatomic, or any specified mixture 
thereof). The code serves the function of the 
tables and charts of NACA Rep. 1135 for any gas 
species or mixture of species, and significantly 
increases the range of valid temperature appli- 
cation due to its thermally perfect analysis. 

Nomenclature 


a speed of sound 

A cross-sectional area of stream tube 

or channel 

Aj coefficients of polynomial curve fit 

for c p /R of individual species; and 
for c p of mixture, equation (3) 

c p specific heat at constant pressure 

c v specific heat at constant volume, 

c p -R 

M Mach number, V / a 

p pressure 

q dynamic pressure, ~ ) - 

R specific gas constant 

T temperature 

u velocity component normal to the 

shock wave 

V flow velocity 

w velocity component tangential to the 

shock wave 

Y mass fraction 

z lateral coordinate to upstream flow, 

measured from wedge leading edge 

y ratio of specific heats, c p / c v 

8 flow deflection angle; 

2-D wedge half-angle 

0 molecular vibrational energy con- 

stant, references 1 and 4 

|i characteristic Mach angle 


p 

mass density 

o 

shock wave angle relative to the 
upstream flow direction 

Subscrij 

its: 

1 

upstream flow reference point; 
e.g., upstream of a shock wave 

2 

downstream flow reference point; 
e.g., downstream of a shock wave 

i 

component gas species of mixture 

j 

coefficient of polynomial curve fit 
for c p 

jmax 

maximum j value 

lim 

limiting conditions for oblique shock 
waves 

max 

maximum value 

mix 

gas mixture 

n 

total number of gas species that 
comprise a gas mixture 

perf 

calorically perfect 

therm 

thermally perfect 

t 

total (stagnation) conditions 

Abbreviations: 

l-D 

one-dimensional 

2-D 

two-dimensi onal 

CFD 

computational fluid dynamics 

CPG 

calorically perfect gas 

cpu 

(computer) central processing unit 

GASP 

General Aerodynamic Simulation 
Program 

NACA 

National Advisory Committee for 
Aeronautics 

NASP 

National Aero-Space Plane 

TPG 

thermally perfect gas 
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Derivation of Oblique Shock 
Relations 

Polynomial Curve Fit for c p 

The selection of a suitable curve fit function 
for c p is the starting point for the development of 
thermally perfect compressible flow relations. 
The form chosen for the TPG code was the eight- 
term, fifth-order polynomial expression given 
below, in which the value of c p has been nondi- 
mensionalized by the specific gas constant. 

C R = ' 4 l(~2) + A 2^) + ^3 + 'Mn + A 5 (T 2 ) (1) 

+ A 6 (r 3 ) + /i 7 (r 4 ) + /i 8 (r 5 ) 


was selected in the current work because of its 
ease of implementation, and the wealth of 
already available data from references 7, 8, and 
9. The only requirements are that closed-form 
solutions to both jc p dT and j(c p /T)dT must be 
known. These requirements are explained in 
more detail in reference 2. 

Mixture Properties 

The TPG code can be used to compute the 
thermally perfect gas properties for not only 
individual gas species but also for mixtures of 
individual gas species (e.g., air). This capability 
is achieved within the code by calculating the 
variation of the heat capacity for the specified 
gas mixture. The mixture c p is 


For a given temperature range, A 2 to A 8 are the 
coefficients of the curve fit. This is the func- 
tional form for c p of NASP TM-1107 (ref. 7), 
which provides values ofA ; to A 8 for 15 different 
gas species. NASA SP-3001 (ref. 8) provides sim- 
ilar c p curve fit coefficient data (A 3 through Ay) 
for over 200 individual gas species for a five- 
term, fourth-order polynomial fit equation. 
NASA TP 3287 (ref. 9) provides additional c p 
curve fit coefficient data for 50 reference ele- 
ments for both a five- term (A 3 -A 7 ) and a seven- 
term (Aj-Ay) fourth-order polynomial fit equa- 
tion. Equation (1) is valid for each of these curve 
fit data sets, with appropriate coefficients set to 
zero. These definitions are typically multi-seg- 
mented curve fits defined over contiguous tem- 
perature ranges so that the c p definition is 
continuous. Thus, for a given species there are 
typically two or more sets of coefficients, each 
valid over a given subset of the entire tempera- 
ture range. 

The polynomial curve fit definitions pro- 
vided in refs. 7, 8, and 9 are defined to be valid 
over specific temperature ranges. Beyond those 
ranges the accuracy of the data fit may degrade 
rapidly. All algorithms based upon these data 
fits should note any calculations at temperatures 
outside of the stated temperature ranges. 

Different algebraic expressions for CpFR 
could be exchanged for that of equation (1) 
within the TPG code. The form of equation (1) 


n 



i = 1 


( 2 ) 


where F; is the mass fraction of the z th gas spe- 
cies. The value of c p is determined from equa- 
tion (1) for each component species. This is 
actually implemented in the code by combining 
like terms of the polynomials for each individual 
gas species; that is 



(3) 


n 

A Jm , x = Z W>, , forj=l,jmax (4) 

/ = 1 

With a known curve fit expression for c p of 
the gas mixture, the value of y for a given tem- 
perature can be directly computed, as can mix- 
ture properties of gas constant and molecular 
weight. See references 2 and 3 for additional 
details regarding the derivation of the 1-D isen- 
tropic relations and the normal shock relations 
from these properties. 

Oblique Shock Relations 

Figure 1 illustrates the geometric relation- 
ships associated with a shock wave at an arbi- 
trary shock angle a relative to the upstream flow 
direction. These relationships can be expressed 
as 
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tano = — (5) 

w 

and 

tan (a - 8) = — (6) 

w 

where Uj and u 2 are the upstream and down- 
stream velocity components normal to the shock, 
and w=wj=w 2 is the tangential velocity compo- 
nent which is conserved across the shock. Since 
only the normal component of velocity changes 
across the shock, the flow is turned through an 
angle 8. Solving equations (5) and (6) for w , a 
single equation can be written: 


tan(o-5) _ u i 

(7) 

tana u x 

The normal 

velocity components can be 

expressed as 


u | = Vjsina 

(8) 

and 


u 2 - V 2 sin(o-8) 

(9) 


Using the ideal gas law ( p = pRT ), equation (13) 
can be written in terms of only temperature and 
normal velocity as 


RT { 
u | 


+ M i 


RT ~ 


(14) 


Substituting for Uj and u 2 with equations (8) 
and (9), equation (14) becomes 


RT } 

V j sina 


+ V j sina = 


RT 2 

V 2 sin (a - 5) 


+ 


V 2 sin(a - 6) , 


(15) 


a function of shock angle o, deflection angle 8, 
the upstream and downstream velocities Vj and 
V 2y and the upstream and downstream static 
temperatures Tj and T 2 . However, as derived in 
reference 2, the velocity may be expressed as a 
function of temperature 


V 2 = 2 ('cdT (16) 

at any point in the flow field, whether upstream 
or downstream of the shock. Thus, the depen- 
dencies on velocity in equations (10) and (15) can 
actually be expressed as dependencies on tem- 
perature. 


Thus, equation (7) can be rewritten as 

tan(p-S) = V 2 sin (q- 5 ) (10) 

tana Vasina * 

a function of shock angle a, deflection angle 8, 
and the upstream and downstream velocities Vj 
and V 2 . 

The continuity and momentum equations 
for 1-D flow across a shock wave in a shock-fixed 
coordinate system (a stationary shock) are given 
as follows 


Assuming all upstream flow conditions are 
known (state 1), and substituting equation (16) 
for V j and V 2 , equations (10) and (15) represent 
a set of two nonlinear equations for three 
unknowns: a, 8, and T 2 . Given any one of these 
variables, the other two may be found numeri- 
cally by means of Newton iteration. (Appendix A 
provides additional details regarding implemen- 
tation of this Newton algorithm within the TPG 
code.) If the shock angle a is known, the tangen- 
tial velocity component is 

w = w | = w 2 = Vjcosa (17) 


Pl"l = P2 M 2 


(11) Thus, the downstream normal velocity u 2 can be 
defined as 


P\ +Pl"] = P2 + P2 W 2 

Dividing the momentum equation by the 
continuity equation gives 

P\ Pi 

h Wi — h U ? 

p,«l p 2 w 2 


( 12 ) 


2 -.2 2 
u, = VS - w 




cdT-w 


(18) 


a function of T 2 only, instead of both T 2 and 8 as 
in equation (9). Thus, instead of equation (15), 
equation (14) reduces to a function of T 2 only 
and can be solved independently of 
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equation (10). For the case of known deflection 
angle 8, equations (10) and (15) must be solved 
simultaneously for T 2 and a. 

Once the primary variables of temperature 
T 2i a, and 8 are known, all other flow quantities 
may be computed. 


V - Pl 

12 C n ~ R 

(19) 

Pl 


a 2 = Y7RT2 

(20) 

M 2 — ^ 7 / ^2 

(21) 


where V 2 is known from equation (16). 

Pressure and density downstream of the 
shock wave can be calculated by the normal 
shock wave relations derived in reference 2, not- 
ing that all velocity terms (including Mach num- 
ber) are defined normal to the shock. That is, 


“ = 1 + y, (A/, sina)^f 


(22) 

P\ v 

«V 



Pi "2 


(23) 


The values of total pressure loss Pt^Ptj and the 
ratio of upstream -static to downstream -total 
pressure p 2 lpt,2 can then be calculated from 
combinations of static and total pressure ratios. 
See reference 2 for further detail. 

Figure 2 illustrates the relationships 
between M h a, and 8 for air at a sample T t . The 
figure shows that, for a given M 2 and 8, two solu- 
tions exist for a, corresponding to what are com- 
monly called ‘weak’ and ‘strong* shock waves. A 
maximum flow deflection angle 8 majc exists for 
any M 2 , defined as the point at which the weak 
and strong shock solutions are identical. Strong 
shock waves always result in subsonic down- 
stream Mach numbers and have shock wave 
angles approaching 90°. Weak shock waves are 
those for which the corresponding wave angle c 
is less than that at which 8 max occurs, and usu- 
ally result in supersonic downstream Mach num- 
bers, except as a approaches o(8 max ). 

Oblique shock waves can only exist at Mach 


numbers above a limiting case, i.e., M 2 >M 2 i im . 
For a given a, the limiting case occurs at p=o, 
where the Mach angle p is defined by 

sinU = tt- (24) 

M \ 

At these conditions the shock strength and 8 
reduce to zero. Obviously, p is only defined for 
M 2 greater than, or equal to, 1, and thus, oblique 
shocks can not exist at subsonic upstream condi- 
tions. For sonic upstream conditions 
M ]—M j iim—i.O, a lim =90° and 6] im =0°. 

For a given 8, the limiting case M 2 j im corre- 
sponds to the Mach number for which 8=8 ma:r 
From a physical perspective, this limit is the 
minimum Mach number at which an attached 
shock wave can occur at the leading edge of a 
wedge. Below M 2 i im the flow results in a 
detached strong shock wave standing a finite 
distance upstream of the wedge leading edge, 
and the flow approaching the point of turning is 
subsonic, not supersonic. 

Oblique Shock Code Description 

An interactive FORTRAN computer code, 
herein referred to as the TPG code, has been 
written based on the equations described herein 
and those described in references 2 and 3. This 
code delivers a complete table of results within 
seconds when run on a computer workstation or 
personal computer. The purpose and primary 
output of the code is the creation of tables of 
compressible flow properties for a thermally per- 
fect gas or mixture of gases (styled after those 
found in NACA Rep. 1135). In addition to the 
isentropic flow and normal shock wave proper- 
ties of the previous version of TPG, the latest 
release of the TPG code (version 3.1) computes 
properties across oblique shock waves, given 
either the shock wave angle a or the flow deflec- 
tion angle 8. Tabular entries may be based upon 
constant increments of Mach number, or con- 
stant decrements of static temperature from the 
total temperature. Such tables are presented in 
NACA Rep. 1135 only for isentropic flow and 
normal shock waves, and only for calorically per- 
fect air (^=1.40) at constant increments of Mach 
number. Only charts (no tables) are presented 
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therein for oblique shock wave properties. As in 
previous versions (see references 2 and 3 ), the 
utility of the TPG code is its capability to gener- 
ate tables of compressible flow properties of any 
gas, or mixture of gases, for any total conditions 
over any specified range of Mach numbers or 
static temperatures ( T<T t , of course). Thus, the 
code recognizes that the properties of thermally 
perfect gases vary with both total temperature 
T t and static temperature T, rather than with 
only the ratio of T/T t . 

In this section, the newest segments of the 
TPG code, Version 3.1, are described, code execu- 
tion procedures are outlined, and examples are 
provided of user-supplied inputs. Various types 
of output are also described. 

The FORTRAN code consists of a Main pro- 
gram, three include files, and 18 subroutines, 
plus a Block Data initialization routine. The 
subroutines are grouped according to purpose 
within nine files and are called as needed from 
various locations within the Main program as 
well as from other subroutines. Version 3 was 
modularized to a greater extent than Version 2 
to simplify and shorten the Main program, and 
to facilitate the internal documentation of the 
source code. Table 1 summarizes the files and 
subroutines which comprise the TPG code. A 
Unix Makefile has been developed for compila- 
tion of the code and creation of an executable. 
Defaults are provided for all inputs and a com- 
plete execution requires, at a minimum, answer- 
ing each interactive prompt with a comma, 
followed by a “carriage return ” Such a minimal 
execution generates a single output file contain- 
ing a table of the basic isentropic flow properties 
for the standard composition of air. More exten- 
sive tables (such as for oblique shock waves), 
tables for other gas mixtures, or files for post- 
processing require specific inputs, which will be 
described in the following subsection. Some 
computer platforms require all character vari- 
able inputs to be delineated by quotes, while 
other computers only require quotes for strings 
which include blank spaces. All platforms on 
which the TPG code has been tested to-date do 
accept quotes. Therefore, for simplicity within 
this report, all character variable inputs are 


Table 1: TPG Version 3.1 files and 
subroutines 


Files 

Subroutines 

TPG.f 

Main program 

Version.h 

Include code version identifier 

params.h 

Include parameters 

dimens. h 

Include common blocks 

initd.f 

Subroutines initd and 
Block Data 

mixture. f 

Subroutine mixture 

trangej.f 

Subroutine trangej 

floprop.f 

Subroutines isntrop and 
shokwav 

cpsubs.f 

Subroutines cpeval , cpNtgra , 
and cptNtgr 

ntgrat.f 

Subroutines intCp and intCpT 

citer.f 

Subroutines Cosigm and 
Codelt 

titer.f 

Subroutines t liter, Todelt , 
Tosigm, sta2lim , and parabt 

postout.f 

Subroutine postout 


shown within quotes. This is the recommended 
procedure for computer platform independence. 

The Main Program 

The Main program handles the interactive 
input and tabular output, and serves as a driver 
for the subroutines which perform the actual 
computation of the mixture properties and flow 
variables. The following paragraphs describe 
the interactive inputs for using the code. Later 
subsections describe the subroutines called from 
Main (see table 1), and provide some details on 
their algorithms and coding. 

The first part of the Main program contains 
the interactive inquiries and responses which 
specify the gas components, the desired total 
temperature, and the desired outputs and out- 
put formats. A typical example of code execution 
is shown in Table 2 in which the user responses 
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are italicized (and numbered in the left margin). 
A response of a comma indicates acceptance of a 
program default. Each user response is followed 
by a carriage return (shown as -i). 

The first response required [1] is an identi- 
fier character string (an ID), or symbol name, for 
inclusion within output files. The ID can be a 
maximum of 70 characters, entered on a single 
line. A default ID (for air) may be specified by a 
“comma” or the word ‘DEFAULT (all upper case 
or all lower case characters), in which case the 
next three questions, [2] -[4], are omitted and 
defaults are assumed. A particular database file 
name may be entered as response [2 7, the pre- 
scribed file containing the necessary thermo- 
chemical data for the gases of interest. The file 
name can be 80 characters, maximum, including 
file path information. The contents of such a 
database file is described in detail in reference 2, 
and is summarized in a later subsection. The 
code contains the database necessary to define a 
four species mixture of air; this default is chosen 
by entering a comma. The gas mixture defini- 
tion is completed by responses [3] and [4] in 
which the number of individual species, and the 
corresponding mass fractions, are specified. 
Defaults are shown as part of the program que- 
ries, and may always be selected; they initially 
correspond to the standard composition of air. 
Regardless of the number of fractions, the code 
checks to see if the sum of the mass fractions is 
equal to unity, and if the database contains suffi- 
cient information for the requested case. If the 
mass fractions do not sum to unity (within a 
small tolerance), the code issues an error mes- 
sage and stops execution. If the sum is within 
unity plus-or-minus the tolerance, the last spe- 
cies fraction is recomputed such that the total is 
precisely unity and a warning message is output 
describing the change. If the database contains 
insufficient information an error message is out- 
put and execution is terminated. 

Responses [1] and [2] are sensitive to the 
case of the alphabetic characters entered. That 
is, upper and lowercase alphabetic characters 
have unique meanings. However, all succeeding 
inputs are insensitive to case, and either upper 
or lowercase responses have the same meaning 


within the code. 

The next inputs define the composition and 
format of the table(s). The total temperature, 
response [ 5 ], specifies the upper temperature 
limit (for which the Mach number is zero). Input 
[6] specifies whether the tables are to be in 
terms of static temperature (T) or Mach (M) 
number increments. Successive entries in the 
table(s) may be given at Mach numbers from an 
initial Mach up to an upper Mach number limit, 
for constant step sizes of Mach, response [7]. 
Alternatively, if the constant static temperature 
option is specified in response [6], the table 
entries are of static temperatures from an initial 
temperature down to the low temperature limit, 
for constant step sizes of static temperature, also 
entered as response [7]. The initial values are 
the total conditions by default. These inputs are 
separated by commas, and any, or all, may be 
defaulted simply by entering a comma alone. 
(See responses [6] and [ 7 7, for example.) 

Response [8] specifies which properties are 
to be included in the table(s). Isentropic flow 
properties are always included in the table(s), 
and constitute the default for response [8]. The 
basic table can be expanded to include properties 
across shock waves through selection of menu 
options 1, 2 or 3. Oblique shock waves can be 
specified through either a specific shock wave 
angle a, or as the shock wave required to turn 
the flow through a given angle 5. Normal shocks 
are a special case of the oblique shock wave with 
o=90°. For menu options 2 or 3, response [9] is 
required to specify the appropriate angle in 
degrees, a or 5, respectively. Response [10] 
allows the user to choose the desired shock wave 
type, either ‘weak’ or ‘strong 5 , for the given 6 of 
menu option 3, response [8], and is omitted for 
all other options. Similarly, responses [9] and 
[10] are both omitted for menu options 0 and 1. 

The relative positions of the tables of isen- 
tropic and shock wave properties are determined 
by response [11] regarding the number of tables 
(either 1 or 2). The shock wave properties and 
the isentropic properties may be output in two 
separate tables cross-referenced by Mach num- 
ber and static temperature, or in a single, wide- 
format table (156 characters across). The sepa- 
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Table 2: Sample TPG v3.1 interactive inputs 


I Thermally Perfect Gas Properties Code I 
I TPG, Version 3.1 .0 I 

I I 

I NASA Langley Research Center I 


Enter the ID ('character string’) for the gas mixture: 

[default = 4-species Air] 

{enter a comma (,) to accept any default} 

[1] ‘Sample Combustion Products’ . I 

Enter the name of the Database File for the Species properties: 

[2] ‘CombProd.data’ J 

Enter the number of species composing the mixture: [default = 4] 

[3] 6 J 

Enter the mass fraction, Y, of all 6 component species 
in the order of the Database File (sum=1 .0): 

[defaults = 0.7553, 0.2314, 0.0129, 0.0004] 

[defaults = 0.0000, 0.0000, 

[4] .1, .6, .2, .0303, .04 J 

Enter the Total Temperature (in degrees K): [default = 400.] 

[5] 1500. J 

Output the table entries by incremental Temp (K) 
or incremental Mach number? (T/M) 

[default = M]: 

[61 , J 

Enter the Initial Mach number, the Mach number Increment, 
and the Upper Mach number Limit: 

[defaults = 0., 0.1, 5.0] 

[7] ,.2,, J 

Enter the appropriate Index from the following Menu: 

0 : Isentropic Flow Properties only, - — > [DEFAULT] 

or additionally include 

1 : Normal Shock Wave Properties, 

2 : Oblique Shock Wave Properties, Shock angle given, 

3 : Oblique Shock Wave Properties, Deflection given, 

[81 3 J 

Enter the Flow Deflection Angle in degrees : [default = 10.000 degrees] 

[91 , J 

Weak or Strong Shock solution: w/s [default = Weak] 

[101 , J 

Enter the number of Tables to be output, 1 or 2 : [default = 2] 

[ 11 ] 1 J 

Enter the Tabular output location: screen=6, file=9 [default = 9 : TPG. out] 

[ 12 ] , J 

Do you want to output a TECPLOT file [TPGpost.dat] ? (y/n): [default = N] 

[131 y J 

Continue calculation to a succeeding state? (y/n) 

[default = N]: 

[14] , J 

STOP: NORMAL 
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rate (second) table of shock wave properties 
always begins with the lowest Mach number for 
which a valid solution of the equations exists, 
that is, Succeeding entries are for Mach 

numbers greater than . In the single table 

format, explanatory text is output under the 
shock property headings when shock property 
data is not valid. 

The table(s) are written to an output file, 
TPG.out, by default since they can be rather 
extensive. The output file can then be viewed 
on-line or printed via an appropriate system 
command. However, if the user desires, the 
tables can be written to standard output (typi- 
cally the computer screen), as specified in 
response [12] to the question of tabular output 
location. 

Response [13] determines the output of an 
additional file for data plotting and postprocess- 
ing. Because of extensive usage at NASA Lan- 
gley Research Center, the Tecplot 1 (ref. 10) code 
format has been chosen for this file. The basic 
postprocessor file TPGpost.dat contains the 
same data as that in the TPG.out tabular out- 
put. The format of this ASCII file, described in 
reference 10 and in appendix B of reference 2, 
could be easily modified to satisfy the input 
requirements of many different graphics post- 
processors. One, two, or three zones of informa- 
tion for all variables and all table entries are 
included. The first zone includes all subsonic 
conditions, and always contains, at a minimum, 
total (stagnation) and sonic entries. The last 
zone contains all supersonic conditions for which 
the desired shock wave exists. For oblique shock 
waves with M lf n m greater than sonic, an inter- 
mediate zone is written containing conditions 
between, and including, sonic and the first shock 
wave. The second and third zones are omitted if 
no supersonic conditions are present. The shock 
wave properties have no physical meaning at 
subsonic conditions (similarly with supersonic 
conditions below the limiting Mach number) and 
are defaulted to a large negative number as a 
warning flag. 

1. Tecplot: trademark of Amtec Engineering, Inc., 
Bellevue, WA 98009-3633. 


Response [14] allows the user to compute a 
succession (or train) of shock waves where the 
upstream conditions of each succeeding shock 
wave is based upon the downstream conditions 
of the preceding one. The default response is to 
end the current set of calculations. An affirma- 
tive response [14] returns the user to the menu 
of options, waiting for a new response [8]. Thus, 
the mixture ID, composition, and total tempera- 
ture are retained, and the table entries are 
determined by the downstream conditions of the 
preceding calculation. 

Following interactive inputs [1]-[13[, the 
database file, if one is requested, is opened, read, 
and closed. Properties of the gas mixture are 
computed and the coefficients of the c p relation- 
ship for the mixture are defined (equations (3) 
and (4)). The TPG.out file is opened and header 
information is written to it stating the gas mix- 
ture specifications and the total temperature. 
The coefficients of the mixture c p relationship 
are also listed (equation (4)). A sample header is 
shown in table 3. 

Through calls from the Main program to 
various subroutines (see table 1) the data arrays 
are initialized, stagnation and sonic conditions 
are computed, and the isentropic flow properties 
are computed for each Mach number or static 
temperature entry in the table, as requested by 
the user in responses [6] and [7], These proper- 
ties are then output to TPG.out unless shock 
properties are requested in conjunction with the 
single table format. If shock wave properties are 
requested (response [8]), the code next deter- 
mines the limiting Mach number and tempera- 
ture conditions for existence of the shock wave. 
Beginning with the first entry at a greater Mach 
number (lower static temperature), the shock 
wave properties are computed to complete the 
data set. Tabular output to TPG.out is com- 
pleted and the file is closed. The same informa- 
tion is then output in postprocessor format to the 
TPGpost.dat file. 

The default “AT (“no”) for response [14] ter- 
minates the calculation. Code execution is com- 
pleted by closing all files and recording warning 
summaries of the times the valid temperature 
ranges were exceeded during the calculations. 
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These warnings, if any, are related to the data- 
base information and are written to an addi- 
tional file called Extrap.WARN for later 
reference by the user. 

Include (or Header) Files 

Three files contain statements that are 
incorporated into the Main program and various 
subroutines via FORTRAN include statements at 
program compilation. The Version.h file con- 
tains a character variable which specifies the 
current version number of the computer code. 
This information is output at the beginning of 
each interactive execution and to the first lines 
of the TPG.out file. 

The params.h and dimens.h files contain 
information which determines the amount of 
computer memory required to execute the TPG 
code, and, thus, the size of the case which may be 
specified. The dimensions of the primary arrays 
within the code as defined below are included in 
a FORTRAN parameter statement in params.h 
and have the following values as released in Ver- 
sion 3.1: 

mtrm (=8): number of polynomial terms 

which define the relationship 
for c p /R y 

mspc (=25): maximum number of gas spe- 

cies allowed, 

mcvf(= 3): maximum number of temper- 

ature ranges per species 
allowed and corresponding 
sets of polynomial coeffi- 
cients, 

mtncr (=10000): maximum number of com- 
puted static temperatures. 

All variables assigned to the preceding 
parameters, except mtrm , may be changed to 
suit particular user requirements. A change in 
mtrm may require coding changes in the subrou- 
tines which perform specific calculations using 
the c p (T)/R relationships. (See the section 
“Polynomial Summation and Integration Sub- 
routines” in reference 2 for details on such 
changes.) The dimens.h file specifies common 


blocks containing all of the major arrays that are 
dependent on the parameters described above. 

Mixture and Basic Flow Properties 
Subroutines 

Subroutine mixture combines the individual 
gas species data using the user-specified mass 
fractions to compute corresponding mixture 
properties. The database file is opened and read, 
if needed. The mixture molecular weight and 
specific gas constant are computed and the mix- 
ture polynomial curve fit coefficients are com- 
puted according to equation (4). The valid 
temperature limits of this mixture curve fit are 
determined from the individual species tempera- 
ture limits. The lower temperature limit of the 
mixture is chosen as the greatest minimum tem- 
perature of all species with nonzero mass frac- 
tion. Similarly, the upper temperature limit of 
the mixture is chosen as the smallest maximum 
temperature of all species with nonzero mass 
fraction. 

The isntrop subroutine computes the basic 
properties of isentropic expansion from total 
(stagnation) conditions for a given static temper- 
ature. The shokwav subroutine computes the 
conditions downstream of a shock wave, relative 
to the upstream conditions from known values of 
the shock wave angle a, the flow deflection angle 
5, the downstream static temperature T 2 , and 
the downstream velocity V 2 . These four funda- 
mental properties are either given as inputs or 
computed in the oblique shock wave subroutines 
( see the following subsection ). 

For any known Mach number M the corre- 
sponding static temperature T can be deter- 
mined from the nonlinear relationship 

2\'c n dT = M 2 yRT (25) 

J J P 

in which y is also a function of T. Subroutine 
tliter performs this function within the TPG 
code. See reference 2 for additional details 
regarding this subroutine and algorithm 
(unchanged from reference 2). 
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Oblique Shock Wave Subroutines 

As stated previously, equations (10) and (15) 
are nonlinear functions of three variables (a, 8, 
and T 2 ) which define the relations across an 
arbitrary oblique shock wave. These equations 
are solved within the TPG code through Newton 
iteration algorithms, given one of the two angles 
a or 5. Subroutine Todelt solves for 8 and T 2 
given a, and Tosigm solves for c and T 2 given 8. 
The downstream velocity component u 2 (and, 
thus, V 2 ) is also obtained in the course of these 
iterations. Each of these iterations is initialized 
by the solution of the calorically perfect shock 
wave equations (see reference 1) using the local 
ratio of specific heats y. Subroutines Codelt and 
Cosigm compute such calorically perfect initial 
estimates. 

Subroutines Todelt and Tosigm only con- 
verge to reasonable answers if the equations are 
physically valid for the given upstream condi- 
tions and desired angles. Thus, the limiting 
static temperature and corresponding Mach 
number Af; /j m must be determined prior to iter- 
ating for specific conditions; subroutine sta2lim 
(short for “station 2 limits”) determines these 
limiting conditions. For a given o, the limiting 
conditions are found at the Mach angle, where 
the shock strength and 8 reduce to zero. For a 
given 8, a two-dimensional iteration procedure 
determines the conditions at which the weak 
and strong shock wave solutions are identical, 
that is, the Mach number at which 5 majc (M)=8. 
The 8 max for a given Mach number is determined 
in subroutine parabt by assuming a parabolic 
form of 8 as a function of a for a given Tp The 
quadratic coefficients of the parabola are 
adjusted until a good fit of the curve apex is 
found, and the resulting quadratic equation is 
solved analytically for the extremum. 

Postprocessor Subroutine 

The postout subroutine in postout.f simply 
opens the TPGpost.dat file, writes the tabular 
data in Tecplot ASCII format, and closes the file 
upon completion. The FORTRAN open state- 
ment specifies a status of “unknown” to allow the 
file to be re-opened for additional output. 


Other Subroutines 

Many of the subroutines in Version 3.1 of 
the TPG code are unchanged since the publica- 
tion of reference 2. Specifically, these routines 
include those which evaluate c p , the integral of 
c p , and the integral of <yT, and are found in files 
ntgrat.f and cpsubs.f. The trangej subroutine, 
which finds appropriate temperature ranges and 
the corresponding polynomial curve fit for a 
given T 9 is also unchanged; it also records 
instances in which the curve fit temperature lim- 
its are exceeded and writes warning messages to 
output files. Finally, the initd and Block Data 
routines are basically unchanged from reference 
2, except that a few isolated variables have been 
added or deleted as the code has evolved. See 
reference 2 for additional details regarding any 
of these routines. 

Database Files 

The thermochemical data required by the 
TPG code for a given mixture of gases is defined 
in a database file to be read at execution time. A 
database for the standard composition of air is 
contained within the code and may be accessed 
as the default. However, for mixtures of other 
gases, or for more complete models of air, a sepa- 
rate database file may be provided. The format 
is simple and modular, grouped by species with a 
two-line header, and additional databases are 
easily constructed. A sample two-species data- 
base is shown in Table 4, and the format is 
described in some detail in references 2 and 3. 

Sample Tabular Output 

Table 5 shows a sample of the tabular out- 
put in the single table format for a shock angle of 
30° in air at 7^=400 K. Following the header 
information, described previously, are columns of 
data for the isentropic flow properties and the 
properties across the oblique shock wave. The 
column headings, in order from left to right, 
refer to the isentropic flow properties of Mach 
number M j (upstream of any shock wave), static 
temperature Tp ratio of specific heats y, static- 
to-total pressure ratio Pj!p t , static-to-total den- 
sity ratio p/ / p£, static-to-total temperature ratio 
Tj!T t , dynamic-to-total pressure ratio q !p t , area 
ratio for expansion from sonic flow A /A*, and 
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local velocity divided by the sonic speed-of-sound 
V/a *. The oblique shock properties continue as 
shock wave angle a, flow deflection angle 5, Mach 
number downstream of the shock wave M 2) 
static pressure ratio across the shock p 2 lpi > den- 
sity ratio across the shock p 2 /p^, static tempera- 
ture ratio across the shock T^/T;, and the ratio 
of total pressures across the shock Pt t 2^Pt,i- For 
a given shock angle, the minimum flow deflec- 
tion is 0°, corresponding to a shock wave of zero 
strength, i.e., a Mach wave. Thus, the minimum 
Mach number of 2.0 for a 30° shock is deter- 
mined by the equation for the Mach angle, equa- 
tion (24). For all non-subsonic Mach numbers 
the code outputs an informative message stating 
the minimum o(=p) for that Mach number. At 
subsonic Mach numbers there is obviously no 
solution to the shock relations, as noted in the 
output. These informative messages do not 
appear in the two-table format; the shock wave 
table simply begins with the minimum Mach 
number Mj n m for which a solution exists. 

Tables 6(a) and (b) illustrate the TPG tabu- 
lar output for a constant flow deflection angle 5. 
Table 6(a) gives the weak shock solutions, and 
6(b) gives the strong shock solutions (i.e., a 
approaches 90° and M 2 <\ ). The column head- 
ings and the subsonic notations are as in table 5. 
However, in this case the limiting Mach number 
is determined by the maximum flow deflection 
angle for which the flow would remain attached 
to the leading edge of a wedge with half-angle 5. 
That is, solutions to equations (10) and (15) 
only exist below some 8 max \ solutions to equa- 
tions (10) and (15) are double valued below 8 max . 
At Mach numbers below the corresponding limit- 
ing Min m , an informative message notes the 
8 max : zero at sonic conditions and increasing 
with Af; to the specified 8. Note that the tables 
always include entries at stagnation, sonic, and 
limiting shock wave M lf u m conditions, even if 
these are not requested by the user. 

Table 7 illustrates the ability of the TPG 
code to calculate flow properties across a 
sequence of shock waves, where each successive 
shock wave is dependent upon the flow condi- 
tions from the preceding one. The simple exam- 
ple consists of Mach 7 free stream flow being 


compressed by two 5° turns. The values of Af 2 in 
the first set of oblique shock wave flow properties 
following the header become the values of M 1 in 
the second set. Version 3.1 treats the shock 
waves independently and generates unique 
table(s) for each. The resultant static pressure 
rise across the two shocks can be obtained by 
multiplying the individual pressure ratios. For 
example, at Af;=7.0: p 2 /pi= 2.235, and 

M 2 =6.104. At M>6. 104: p 2 / P] =2.034 , and 
M 2 - 5.398. The resultant static pressure rise is 
therefore 2.235 x 2.034 = 4.546. The other prop- 
erties of density, temperature, and total pressure 
ratios across both shock waves are calculated in 
a similar fashion. A future version of TPG may 
condense the data into a single table and provide 
cumulative ratios automatically. 

Oblique Shock Code Validation 

The general polynomial curve fit methodol- 
ogy of the TPG code for thermally perfect proper- 
ties was previously validated in references 2 and 
3 . The oblique shock wave capabilities within 
the TPG code were validated in the same man- 
ner as were the basic isentropic flow and normal 
shock wave methods in these references, both in 
the calorically perfect and in the thermally per- 
fect temperature regimes. In the current devel- 
opment, a normal shock wave is simply a special 
case of an oblique shock wave with g=90°; the 
previous validations still hold for this special 
case. 

Calorically Perfect Temperature 
Regime 

The first validation test of the TPG code’s 
oblique shock wave capability was the verifica- 
tion of accuracy in the temperature regime usu- 
ally considered to be calorically perfect, with air 
as the test gas. In this temperature regime, the 
specific heat of air is nearly constant and there- 
fore the TPG code results should be nearly iden- 
tical to results obtained from the calorically 
perfect formulas of NACA Rep. 1135. For these 
test cases, the default data file for standard four- 
species air was used along with the default 
(standard) values for the mass fraction composi- 
tion of air. The code default total temperature of 
400 K was selected, a temperature considered 
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within the calorically perfect temperature 
regime for air. The lower Mach number limit 
was defaulted to stagnation conditions and the 
upper Mach number limit was set to 10, with a 
Mach number increment of 0.1. In the case of 
real air at 7^=400 K, liquefaction would occur 
well before Mach 10; the data is presented to 
that extreme herein merely for comparison with 
the calorically perfect gas tables of NACA Rep. 
1135 (which, incidentally, extend up to Mach 
100). (All cases presented herein were run with 
the TPG code in the Mach number increment 
mode. See reference 2 for examples of the tem- 
perature increment mode.) 

Figures 3(a-c) and 4(a-c) compare TPG cal- 
culations of M 2 , 8, p 2 /Pi, T 2 /T h p 2 /p h andp t2 /p t j 
for air with calorically perfect (y=1.40) solutions 
at shock wave angles of 30° and 50°, respectively. 
Each variable is plotted versus upstream Mach 
number Mj. The differences are negligible as 
shown. Figures 5 and 6 show the same data pre- 
sented as ratios of the thermally perfect values 
to the calorically perfect values. At a=30°, the 
differences are less than 0.23% for all variables 
except total pressure ratio Pt^Pt 1- Similarly, for 
<7=50° the differences are less than 0.25%, except 
for total pressure ratio. Even for this most sensi- 
tive variable, the differences are less than 1%, 
and that value is approached only at the largest 
shock angle and for Mach numbers approaching 
10. As noted in reference 2 for normal shock 
waves, these small differences actually represent 
the small amount of error associated with the 
calorically perfect gas assumption. A slight vari- 
ation with temperature actually exists in the 
heat capacities even at such a low T t) and is the 
cause of the variations seen in figures 5 and 6. 

Thermally Perfect Temperature 
Regime 

The properties for arbitrary oblique shock 
waves were verified for test cases involving stan- 
dard four-species air at total temperatures of 
1500 K and 3000 K, and at shock wave angles of 
30° and 50°. Comparisons were made with 
results obtained via the imperfect gas relations 
of NACA Rep. 1135 (i.e, the 0-equations for a 
diatomic gas) and with selected results from the 
calorically perfect gas assumption. The proper- 


ties of 5, M 2 , p 2 /pi, T 2 /T h p 2 /p/, and Pt^Pt.i are 
plotted versus the upstream Mach number Mj in 
figures 7 and 8. Figures 7(a)-(f) show results for 
7^=1500 K, a temperature well above the limits 
of the calorically perfect assumption for air, and 
figures 8(a)-(f) are for 7^=3000 K, an even more 
extreme case. The two thermally perfect gas 
methods agree extremely well. Cursory exami- 
nation of the plots shows that imperfect gas 
effects increase with both total temperature and 
shock angle. For some flow conditions the differ- 
ences appear to be negligible, but, for other con- 
ditions, the differences can be large. Also, not all 
properties appear to be equally sensitive to 
caloric imperfections. 

The calculation for a total temperature of 
3000 K is presented as an extreme case to illus- 
trate application of the TPG code over an 
extended temperature range under the assump- 
tion of chemically frozen flow, i.e., fixed composi- 
tion. The presentation is in the same spirit as 
that of references 5 and 6 f where data is pre- 
sented to 5000 K and 6000 K, respectively. As 
always, users must keep in mind the applicabil- 
ity of the assumption of chemically frozen flow to 
a particular problem. 

A more descriptive mode for examining the 
differences between the thermally perfect gas 
calculations and calorically perfect results is 
shown in figures 9(a)-(f), In these plots, the 
thermally perfect properties (8, M 2 ,p 2 !p h T 2 /T h 
p 2 /p;, and Pt^Ptj} are divided by the calorically 
perfect properties at identical values of Mj to 
yield a measure of the imperfect gas effects. 
Results are shown for both total temperatures 
and both shock angles, for both the TPG code 
and the 0-equation method of NACA Rep. 1135 
(ref. 1). In this format differences due to caloric 
imperfections are observed of 10% or greater for 
temperature ratio, density ratio, and total pres- 
sure ratio at the higher total temperature and 
greatest shock angle. Caloric imperfection dif- 
ferences on the order of 4-5% are observed for all 
variables except static pressure ratio, even at 
the lower total temperature. Also, the magni- 
tudes of the caloric imperfections do not all vary 
uniformly with Mach number M ; . The effects on 
8 are largest at lower Af;, but are larger for other 
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variables at higher Mj. For static pressure ratio 
the sign of the caloric imperfection effect 
changes as M 1 increases. Thus, the use of the 
calorically perfect relationships could result in 
substantial errors, particularly in the design or 
analysis of flows with a sequence of such shock 
waves. Finally, these results are shown assum- 
ing that Mj is the control variable, i.e., the vari- 
able held fixed for both calorically and thermally 
perfect calculations. For other applications, the 
control variable might be the flow deflection 
angle 8, for example, and the magnitude of the 
caloric imperfection could be even greater. 
(Comparisons relative to the total, state 0, condi- 
tions, rather than state 1, would also yield larger 
effects of caloric imperfection.) 

The differences between the two thermally 
perfect gas methods can be seen to be small (typ- 
ically less than 0.25%). This excellent agree- 
ment verifies the accuracy of the TPG code with 
the derived thermally perfect gas relations for 
oblique shock waves based upon polynomial 
expressions for c p . Although these test cases 
were all for air, which is a primarily diatomic 
gas, the TPG code is also valid for polyatomic 
gases. The utilization of a polynomial curve fit 
for c p makes the TPG code applicable to any 
molecular structure of the gas, not just the 
diatomic structure upon which the ©-equations 
are based. 

Comparisons with CFD 

As an additional validation of the TPG 
code's oblique shock wave capability, results 
were compared to computational fluid dynamics 
Euler solutions (i.e., in viscid flows). The CFD 
solutions were obtained using the General Aero- 
dynamic Simulation Program (GASP), ref. 11 , of 
AeroSoft, Inc. GASP solves the integral form of 
the governing equations, including the full time- 
dependent Reynolds-averaged Navier-Stokes 
equations and various subsets: the Thin-Layer 
Navier-Stokes equations, the Parabolized 
Navier-Stokes equations, and the Euler equa- 
tions. A generalized chemistry model and both 
equilibrium and non-equilibrium thermodynam- 
ics models are included. The discretized equa- 
tions may be solved by space marching or by 
global iteration. The current computations uti- 


lized space marching (i.e., totally supersonic 
flow) for the inviscid Euler equations. Third- 
order upwind inviscid fluxes were calculated 
using Roe’s split flux normal to the flow direction 
with Min-Mod flux limiting. The NASA Lewis 
Research Center equilibrium curve fits for spe- 
cific heat (references 7, 8 and 9), of the same 
form employed within TPG, defined the thermo- 
dynamic characteristics of the gases. A 2-D grid 
of 20 uniformly-spaced computational cells per 
unit length was defined on a nondimensional 
domain (10 units long) above a simple wedge. 
The domain height was chosen so that the 
resulting oblique shock wave would pass 
through the downstream boundary (5 and 8.4 
units for the two test cases, respectively). 

The first CFD test case involved conditions 
representative of a generic high-speed inlet: 
thermally perfect air at Mj= 4.8, Tj=l610 K and 
/?l=101325 N/m 2 (1 atmosphere) flowing over a 
10° wedge. Figure 10 shows contours of pressure 
ratio pip i for the GASP Euler solution with the 
shock wave computed by the TPG code 
(a= 19.415°) overlaid as a thick dashed line. The 
magnified view (see the inset) shows that the 
agreement in shock position is excellent, consid- 
ering that GASP captures the shock over a num- 
ber of cells and the TPG solution is a precise 
point-solution, that is, the shock wave is a true 
discontinuity. More detailed comparisons of M, 
pip TIT ly and p/p; are shown in the line plots of 
figure ll(a)-(d). The data points shown were 
interpolated from the GASP solution field (e.g., 
figure 10) along a horizontal line at z- 2.6 units 
above the wedge leading edge. The TPG solution 
is shown as a discontinuous solid line with the 
shock position calculated geometrically at that 
same lateral location of 2.6 units. Again, the 
agreement is excellent outside the region of 
shock capturing. 

The ©-equations of NACA Rep. 1135 are 
valid only for a single-species diatomic gas. 
Therefore, a second comparison with GASP was 
computed at arbitrary conditions designed to 
highlight the flexibility of the TPG code. The gas 
chosen was a two-species mixture of 20% steam 
and 80% C0 2 by mass, neither species being 
diatomic. Upstream conditions were set at 
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Mj= 3.5, 7>1200 K and /> ; =101325 N/m 2 . A 
rather severe wedge angle of 25° was selected 
and the 2-D GASP grid was constructed as 
before: 20 uniformly-spaced computational cells 
per unit length defined on a nondimensional 
domain (10 by 8.4 units to contain the expected 
shock wave). GASP required 282.4 cpu seconds 
to compute the solution on a Cray 2 Y-MP com- 
puter, not including grid generation and GASP 
input setup time (about 2 days). The TPG solu- 
tion was practically instantaneous, computed on 
a Sun 3 Sparcstation 20 workstation. The result- 
ing solutions are shown in figure 12, again with 
the TPG solution (a=37.3037°) overlaid as a 
thick dashed line on the GASP pressure ratio 
contours. The magnified view in the inset again 
confirms the excellent agreement of the shock 
positions. The line plots of M, p/p h T!T h and 
p/p; shown in figure 13(a)-(d) confirm the accu- 
racy of the TPG code for arbitrary mixtures of 
polyatomic gases. Since the TPG result is a pre- 
cise point-solution, the overshoots seen in 
figure 13, and common to shock-capturing CFD 
codes, are eliminated. Similar to figure ll(a)-(d), 
the data points were interpolated from the 
GASP solution (figure 12) along a horizontal line 
at z- 5.0 units above the wedge leading edge, and 
the TPG solution is shown as a discontinuous 
solid line with the shock position calculated at 
the same lateral location of 5 units. 

Temperature Limits of the TPG 
Code 

The TPG code provides valid results as long 
as its application is within the thermally perfect 
temperature regime for the gas of interest. If 
the code is used outside of the region for which c p 
is a function of temperature only, then the ther- 
mally perfect results produced by this code will 
no longer accurately reflect what actually hap- 
pens in nature. Dissociation will occur above 
some upper temperature limit, and will cause a 
deviation from thermally perfect theory. At the 
opposite extreme, below some lower temperature 
real gas effects will become important (i.e., inter- 
molecular forces will not be negligible). The 

2. Cray Research, Inc., Minneapolis, MN 55402. 

3. Sun Microsystems, Inc., Mountain View, CA 
94043. 


TPG code user must remain aware of these high 
and low temperature boundaries associated with 
the thermally perfect assumption for the partic- 
ular gas(es) under consideration. Note that 
these are boundaries associated with the physi- 
cal properties of the gas(es), and are distinct 
from the upper and lower temperature limits 
associated with the polynomial curve fits utilized 
within the TPG code. These latter curve fit lim- 
its are tracked within the code, and users are 
warned by the TPG code when these limits are 
exceeded. The limits of the thermally perfect 
assumption are the user’s responsibility, and are 
dependent upon the specific gas properties and 
flow conditions. See reference 2 for a more 
detailed discussion of these temperature limits. 

Conclusions 

A set of compressible flow relations describ- 
ing flow properties across oblique shock waves 
has been derived for a thermally perfect, calori- 
cally imperfect gas, and applied within the exist- 
ing thermally perfect gas (TPG) computer code. 
The relations are based upon a value of c p 
expressed as a polynomial function of tempera- 
ture. The code now produces tables of compress- 
ible flow properties of oblique shock waves, as 
well as the properties of normal shock waves and 
basic isentropic flow from the original code, in a 
format similar to the tables for normal shock 
waves found in NACA Rep. 1135. The code 
results were validated in both the calorically 
perfect and the calorically imperfect, thermally 
perfect temperature regimes through compari- 
sons with the theoretical methods of NACA Rep. 
1135, and with a state-of-the-art computational 
fluid dynamics code. The TPG code is applicable 
to any type of gas or mixture of gases and is not 
restricted to only diatomic gases as are the ther- 
mally perfect methods of NACA Rep. 1135. 

The TPG code computes properties of 
oblique shock waves for given shock wave angles, 
or for given flow deflection (wedge) angles, 
including both strong and weak shock waves. 
Normal shock waves are a special case of an 
oblique shock wave. The code also allows calcu- 
lation of a sequence of shocks where each succes- 
sive shock wave is dependent upon the flow 
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conditions from the preceding one. Both tabular 
output and output for plotting and post-process- 
ing are available. Typical computation time is 
on the order of 1-5 seconds for most computer 
workstations or personal computers. The user 
effort required is minimal, and consists of only 
simple answers to 14 interactive questions (or 
less) per table generated. Thermally perfect 
properties of flows through oblique shock waves 
may be computed with less effort than has tradi- 
tionally been expended to compute calorically 
perfect properties. Errors incurred from using 
calorically perfect gas relations, instead of the 
thermally perfect gas equations, in the ther- 
mally perfect temperature regime have been 
shown to approach 10% for some applications. 

The TPG code is available from the NASA 
Langley Software Server (LSS) on the World- 
Wide- Web at the following address: 

http://www.larc.nasa.gov/LSS 
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Appendix A 


and equation (A7) is actually 


Equations for Numerical 
Implementation 

Equations (7) and (14) may be written as 

F,(o,8, T 2 ) = tan (q- 5 . ) _^ = Q ( A1 ) 

J z tana u ] 


and 


G 2 (S, T 2 ) = G 2 (T 2 ) = 0 (A9) 

Once equation (A9) is solved for T 2 , V 2 is known 
from equation (A4), u 2 from equation (A3), and 8 
may be found from equation (A6) or simply from 

cos(a-S) - (A10) 

V 2 


RT , RT 1 

F 2 (a 9 5, T 2 ) - — + Ui - = 0 (A2) 

u \ u 2 

where state 1 conditions T 1 and uj are assumed 
to be known, and u 2 is a function of a, 5, and T* 
as seen from equations (9) and (16), reproduced 
here for completeness. 


In the case of the flow deflection angle 8 
being given, equations (Al) and (A2) are a sys- 
tem of two nonlinear equations for the two 
unknowns a and T 2 , and may be written as 


T 2 ) 


tan(q - 8) 
tana 


V 2 sin(a - 8) 
V ] sin a 


(All) 


u 2 = V 2 sin(a-8) 




2 

2 



(A3) 

(A4) 


and 

F 2 (a, T 2 ) 


RT ] RT 2 

: — + V, sina : =- 

VjSina 1 V / 2 sin(a-8) 

-V 2 sin(a - 8) (A12) 


Thus, there are two nonlinear equations for 
three unknowns. Typically, either a or 8 is 
known, and the equations may be solved numer- 
ically for the remaining angle and T 2 by Newton 
iteration. 

In the case of the shock wave angle a being 
given, the equations may be rewritten as 

F. (a, 8, T 2 ) = G.(8, r 2 )| o , for i=l,2. (A5) 

That is, 

G, (8, T 2 ) = tan(g ~ S) _— = 0 (A6) 

1 2 tano m, 

and 

RT, RT - 

G 2 (8, T 2 ) = i + - u 2 = 0 (A7) 

u \ u 2 

Now u 2 can also be written as 


The solution can be found by application of New- 
ton’s method to a system of equations, similar to 
the application to equation (A9). 

For equation (A9), a single equation of a 
single unknown, Newton’s iterative method can 
be written 


(A7\,)* +] +G* = 0 (A13) 

where the superscripts denote the iteration step. 
That is, the function G 2 and the derivative of G 2 
are evaluated using values at the current, or 
k’th , iterative approximation step in order to 
solve for the k+l’th correction to the approxi- 
mate solution. The differential equation is then 
replaced by the following algebraic equation 



(A7- 2 )* + 1 = 7-* + ‘ -T k 2 = C 2 



(A14) 


2 

u 2 




c p dT - w 2 


(A8) 


Therefore, u 2 is a function of only T 2 if a is given, 


For a system of equations, the functional 
evaluation G 2 * becomes a vector, and the deriva- 
tive becomes a matrix of partial derivatives 
which must be inverted. The correction term 
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A T 2 becomes a vector of the unknowns o and T 2 , 
in the case of equations (All) and (A12). Equa- 
tion (A14) is thus replaced by a system of two 
algebraic equations for the two unknowns. 

Newton’s method requires an initial esti- 
mate of the solution. Within the TPG code, the 
initial approximation is given by the solution of 
the oblique shock wave relations for a calorically 
perfect gas at the local conditions (i.e., Afj, and 
7 ;). NACA Rep. 1135 summarizes these equa- 
tions, as do many compressible flow textbooks. 

Newton’s method applied to equations (A9) 
and (All-12) does not converge for M 1 <M h i im 
where oblique shock waves cannot exist, as dis- 
cussed previously in the section on Oblique 
Shock Relations . For the case of known a, the 
limiting Mj is defined by equation (24), where 
I i=g. For the case of known 5, the limiting Mj is 
defined along the curve labeled S max in figure 2. 
Within the TPG code, an additional nonlinear 
iteration searches for the point on the curve at 
which 8 ma;t =5. At each step of the search, the 
value of b max for a given Mj is found by parabolic 
approximation of the curve apex. 
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Table 3: Sample TPG.out header information 


Thermally Perfect Gas Properties Code 
TPG, Version 3.1.0 

NASA Langley Research Center 


Table of Thermally Perfect Compressible Flow Properties 
Gaseous Mixture: Air: 4-Species Mixture of N2, 02, Argon, and C02 
Database file name: [ Default database of Air Mixture ] 

Species Names are: 

N2 , 02 , Argon, C02 , 

Species Mass Fractions are: 

0.75530 0.23140 0.01290 0.00040 

Species Mole Fractions are: 

0.78092 0.20946 0.00935 0.00026 

Mixture Properties: 

Molecular Weight = 28.9663 

Gas Constant = 2.87035E+02 J/(kg*K) 

Polynomial Coefficients : cp/R = Sum{A(i)*T A (i-2)} 

80.0 < T < 1000.0 degrees K 

-3.731 9821 D+00 5.41 96701 D-01 3.4693586D+00 3.8353798D-04 

-2.76531 30D-06 8.1886142D-09 -7.7368690D-12 2.4348723D-15 

1000.0 < T < 6000.0 degrees K 

2.3769851 D+05 -1 2440527D+03 5. 1 266690D+00 -2.0400644D-04 

6.8321 801 D-08 -1 0553542D-1 1 6.6450071 D-1 6 0.0000000D+00 

Total Temperature = 400.000 K 
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Table 4: Sample thermochemical database file 


# of species #of Temperature ranges : NASP TM 1107 + Mods. 

2 2 
Name: Molecular Wt. 

N2 28.0160 

tmin tmax 

20 . 1000 . 

Cp/R Coefficients: cl (-2) — > cl (5) 

-1.33984200E+01 1 .34280300E+00 3.45742000e+00 5.74727600E-04 

-3.2171 1900E-06 7.50775400E-09 -5.90150500E-12 1 .50979900 E-1 5 

tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2(-2) — > c2(5) 

5.87702841 E+05 -2.23921 563E+03 6.06686971 E+00 -6.13957913E-04 

1 .491 78026E— 07 -1 .923071 30E-1 1 1 .06193594E-15 0.00000000E+00 

Name: Molecular Wt. 

02 32.0000 

tmin tmax 

30. 1000. 

Cp/R Coefficients: cl (-2) — > cl (5) 

3.8851 7500E+01 -2.70630800E+00 3.561 19600E+00 -3.32782400E-04 

-1.18148000E-06 1 .10853500E-08 -1 .49299400E-1 1 5.99553800E-1 5 

tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2(-2) — > c2(5) 

-1 .05642070E+06 2.41 123849E+03 1 .73474238E+00 1.31512292E-03 

-2.299951 51 E-07 2.13144378E-1 1 -7.87498771 E-1 6 0.00000000E+00 
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Table 5: Sample TPG tabular output for constant shock wave angle 


I Thermally Perfect Gas Properties Code j 
t TP'j, Version ?.1.0 

t I 

I NASA Langley Research Center ( 


Table of Thermally Perfect Compressible Flow Properties 
Caseous Mixture: Air: 4 -Species Mixture of Nr. 02. Argon, and C02 

Database file name: [ Default database of Air Mixture j 
Species Names are: 

N2 . (.(? , Argun. CO? . 

Species Mass Fractions are: 

0.75530 0.23140 0.01290 0.00040 

Species hole Fractions are: 

0.78092 0.20946 0.00935 0.00026 

Mixture Properties: 

Molecular Weight * 28.9663 

Cas Constant - 2.870356 + 02 T/Cfccft.) 

Polynomial Coefficients : cp/fi 1 Sui-fACi J*f A (i -2)} 

80.0 T 1000.0 deques K 

-3. 7319821 D+UO 5.4196701 D-01 3. 469 35860+00 3.8353798D-04 

-2. 76571 30D-I3E 8. 1806142D-09 -7. 73686900-1 2 2.43487230-15 

KiOO.O T 6000.0 degrees K 

2. 37698510+05 -1 . 24405 27()+ u? 5. 12686900*00 -2. 04 00644 ['-04 

6.8321801D-08 -1 . 05535420- 11 6.64500710-16 0. 00000000+00 

Total Temperature = 400.000 F 


ISENTR0PIC FLOW 

PROPERTIES ----- 






OBLIQUE SHOCK 

WA 

iT FLOW PROPERTIES (2-downstream, 

1 -upstream) — 

M 

T ( h) 

Gamma 

P/Pt 

RHO/RHOt 

T/Tt 

q/Pt 

A/A ; 

V/a* 

ShoFAngl FluDfl 

M2 

P2/P1 RH02/RH01 

T2/7 1 

PfLVPtl 

0.000 

400.0 

. 3949 

1. OOOE+OO 

1 . 000E+00 

1 . nnnn 

0.0006+00 

infinite 

0.00000 

+ * Subsur 


Net Appl 

cable — 



0.200 

396.9 

. 3951 

9.726E-01 

9. 803E-01 

0. 9922 

2.7146-02 

2. 966E+GG 

0.21784 

* ; Sub sen 


Not App’M 

cable --- 



0.400 

387.7 

.3956 

8. 959E-01 

9.242E-01 

0.9693 

1.000E-01 

1.591E+00 

0.43073 

**• Subsnrnc 

Nc-r Appl 

cable — 



0.600 

373.4 

. 3965 

7. 844E-G1 

6.403E-01 

0. 9335 

1.9726-01 

1 . 1 89E+QQ 

0.63422 

++ Subsonic 

Net Appl 

cable — :■ 



0.800 

355.0 

.3975 

6.564E-01 

7.397E-Q1 

0. 8674 

2.935E-01 

1 . 038E+00 

0.82478 

+ n suhsom c 

Not Appl i cabl e — 



1.000 

333.7 

. 3985 

5. 285E-01 

6 . 335E-01 

0.8343 

3.695E-01 

1 . 000E+00 

1.00000 

*■- Minimum 

toMque Shck 

k Wave Angle - 90.000 



1 . 200 

310.9 

.3994 

4.124E-01 

5. 305E-01 

0,7773 

-4. 1 556-01 

1.031E+00 

1.15866 

+ + Minimum 

ttlique Shock Wave Angle = 56.443 



1 .400 

287.7 

.4001 

3.141E-01 

4. 368E-01 

0.7192 

4.310E-01 

1 . 115E+00 

1.3005? 

1, " t: Minimum 

ttlique SHck 

k Wave Angle - 45.585 



1.600 

264.3 

.4006 

2. 351 E— Cl 

3.551 E-01 

0.6620 

4. 2146-01 

1 .251E+00 

1.42631 

** Minimum 

•b 1 i que Sho< 

k Wave Angle = 38.682 



1.000 

242.9 

.4010 

1.739E-01 

2. 863E-01 

0.6072 

3.946E-01 

1.440E+00 

1.53700 

Minimum 

til ique ShcK 

k Wave Angle ■ 33.749 



2.000 

222.3 

.4012 

1.277E-01 

2.2976*01 

0.5558 

3.577E-01 

1 . 688E+00 

1 .63404 

30.0000 0. 

001 

0 2.000 

1 . 0006+00 1 . OOOE+OO 1 

OOOE+OO 

1 .OOOE+OO 

2.200 

203.3 

.4013 

9.341E-02 

1 .838E-01 

U . 508 3 

3. 168E-01 

2. 005E+00 

1 . 71892 

30.0000 3. 

7160 2.059 

1 . 245E+00 1 . 1 69E+00 1 

. 065E+00 

9.983E-01 

2.400 

185.5 

.4013 

6.833E-02 

1 .4706-01 

0.4648 

2.758E-01 

2.403E*Q0 

1 .79311 

30.0000 6. 

71 

3 2.130 

1.514E+00 1 . 34 IE +00 1 

. 1 28E+00 

9. 928E-01 

2.600 

1 70. 1 

.4013 

5.008E-02 

1 .1786-01 

0.4252 

2.3726-01 

2.895E+00 

1.85799 

30.0000 9. 

14 . 

1 2.208 

1.805E+00 1.5156+00 1 

.1916+00 

9. 794E-01 

2.800 

155.0 

.4013 

3. 684E-D2 

9.46DE-02 

0. 3894 

2.023E-01 

3.4976*00 

1.914 80 

30.0000 11. 

12? 

1 2.28? 

2.120E+00 1.689E+00 1 

. 255E + 00 

9. 582E-01 

3.000 

142.8 

.4013 

2. 722E-02 

7. 624E-02 

0.3571 

1 . 71 76-0 1 

4.2296*00 

1.96463 

30.0000 12. 

?6 C 

6 2. 367 

2.459E+GG 1.861E+00 1 

. 321E+00 

9. 298E-01 

3.200 

111. 2 

.4012 

2. 024E-02 

6. 169E-02 

0.3280 

t. 4526-01 

5.1 12E+G0 

2.00845 

30.0000 14. 

12 ? 

0 2.444 

2. 821E+00 2.031E+00 1 

. 389E+00 

8. 953E-01 

3.400 

120.7 

.4012 

1 .514E-02 

5. 015E-02 

0.3018 

1 . 226E-U1 

6.171E+00 

2.04 707 

30.0000 15. 

27 

4 2.520 

?. 206E+00 2.196E+00 1 

.460E+C0 

8.558E-01 

3.60C 

111.3 

.4012 

1.140E-02 

4.096E-02 

0.2783 

1 .O35E-0I 

7.431 E*00 

2.08122 

30.0000 16. 

23 ' 

1 2.592 

3. 614E+00 2. 357E+00 1 

.533E+00 

8. 126E-01 

3.800 

102.8 

.4012 

8.E43E-03 

3. 362E-02 

0.2571 

8. 743E-03 

8.9246*00 

2. 1 1 149 

30.0000 1". 

0643 2.661 

4.04GE+00 2.514E+00 1 

.6106+00 

7.676E-01 

4.000 

95.2 

.4012 

6.599E-03 

2.773E-02 

0.2380 

7. 397E-02 

1.0686+01 

2.13842 

30.0000 17. 

77 

3 2.727 

4. 501E+00 2 . 664E+00 1 

. 689E+00 

7. 21 IE-01 
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Table 6: Sample TPG tabular outputs for constant flow deflection angle 


(a) Weak shock wave solutions 


l Thermally Perfect Gas Properties Code i 
\ TPG. Version 3.1.0 ! 

I I 

t NASA Langley Research Center I 


Table of Thermally Perfect Compressible Flow Properties 

Caseous Mixture: Air: *4 -Species Mixture of HI'. 02, Argon, and C02 

Database file name: [ Default database of Air Mixture ] 

Species Names are: 

N2 , 02 , Argon. C02 , 

Species Mass Fraction's are: 

0.75530 0.23140 0.01290 0.00040 

Species Mole Fractions are: 

0.78092 0.20948 0.00935 0.00026 

Mixture Properties: 

Molecular Weight = 28.9663 

Gas Constant - 2.87035E+02 J/(kg*K) 

Polynomial Coefficients : cp/k * Sum; A(i )*T' k (i -2)} 

80.0 ' T 1000.0 degrees K 

- 3. 7 3 14821 [>+00 5.41 967010 -CM 3 .4693586P+00 3. 8 35 37960-04 

-2.76531 300 -06 8 . 1 8861428-09 -7. 7368690D-1 2 2.43487230-15 

10ijU.li T 6000. n degrees h 

2. 3769851 [i ♦OS -1 . 2440527D+O2 5.12666900*00 -2. 0400644D-04 

6. 8 321 801 P-GO 1 . 05535471'- 1 1 b. 64500710- 16 0. 00000000+00 

Total Temperature « 400.000 K 


ISENTR 

OPIC FLOW 

PROPERTIES — 






OBLIQUE SHOCK WAVE FLOW PROPERTIES 

2 =downstream , 1 =upstreait 

} 

M 

T!Kf 

Gamma 

D/Pt 

RHO/RHOt 

T/Tt 

q/Pt 

A/ A* 

V/a f 

ShokAngl FloDflec M2 

P2/P1 

RH02/RH01 T2/T1 

Pt2/Pt1 

0.000 

400.0 

. 3949 

1 . O0OE + OQ 

1 .000E+00 

1 . 0000 

0. OOOE+OG 

i nfini te 

0 . 00000 

♦ + Subsonic: Not App 

li cable — 



0.200 

396.9 

.3951 

9. 72RE.-01 

3.003E-01 

0.9922 

2 . 7 HE -02 

2. 966E+00 

0 . 21 784 

* + Subsonic: Not App 

licable — 



0.400 

387.7 

. 3956 

8.959E-0I 

9. 242E-01 

0.9693 

1 . 0Q0E-O1 

1 . 591E+00 

0.43073 

+ + Subsonic: Not Applicable — 



0.600 

373.4 

. 3965 

7.844E-Q1 

8.403E-01 

0.9335 

1 .972E-01 

1.189E+00 

0.63422 

*+ Subsonic: Nut App 

n cable -• 



0.600 

356.13 

.23975 

6.564E-01 

7.397E-01 

0.8874 

2.935E-01 

1.038E+00 

0.82478 

** Subsonic: Not App 

licable --- 



1 . 000 

333.7 

. 3985 

5.265E-01 

6 . 335E-01 

0.8343 

3. 695E-01 

1. 000E+00 

1 . 00000 

*+ Maximum Attached-f 

low Shock be 

flection Angle * 0.000 


1 . 21.10 

310.6 

. 39*34 

4. 124E-01 

5.305E-01 

0.77721 

4. 155E-01 

1 . 1)31 F +00 

1 . 1 5866 

* + Maximum Attached-f 

low shock [reflection Angle - 3.95'! 


1 .400 

287.7 

.4001 

3.141E-01 

4. 368E-01 

0.7192 

4. 310E-01 

1.11 5E+00 

t . 30057 

* + Maximum Attached-f 

low Shock Deflection Angle = 3.444 


1 . 600 

264.6 

. 4006 

2. 351E-01 

3.551E-01 

0.6620 

4.214E-01 

1 . 25 IE +00 

1 .4263: 

+ + Maximum Attached-f 

low Shock Deflection Angle » 14.671 


1 . 800 

242.9 1 

.4010 

1 . 739E-Q1 

2 . 863E-01 

0.6072 

3. 946E-0I 

1 . 440E+00 

1 . 53700 

* + Maximum Attached-f 

low Shock Deflection Angle = 19.201 


1 . 639 

238.7 1 

. 411 1 0 

1 .637E-01 

2. 743E-01 

0.5969 

3.879E-01 

1.483E+00 

1.55704 

64.9121 20.0000 0.920 

3. 073E+OD 

2.M3E + 0U 1.434E+U0 8 

.696E-01 

2.000 


.4012 

1.277E-01 

2.297E-01 

0.5558 

3.577E-01 

1 . 688E+00 

1 .63404 

53.4172 20.0000 1.211 

2. 844 E +00 

2. 042E+00 1 . 393E+00 5 

.329E-01 

2.200 

203.3 1 

.4(31 3 

9. 341 E- 02 

1 .838E-G1 

0.5083 

3. 168E-01 

2. 005E+00 

1.71892 

47.9823 20.0000 1.403 

2. 951 E+ 00 

2.089E+00 1.41 3E + 00 8 

.82 IE -01 

2.400 

185.9 1 

. 40 1 5 

6. 83 3E-02 

I .470E-01 

0.4648 

2. 758E-01 

2.403E+00 

1.793M 

44.3478 20.0000 1.568 

3.118E+00 

2.160E+0G 1.444E+00 6 

.650E-01 

2.600 

170.1 1 

. 40 1 2 

5. 008E-D2 

1 . 178E-01 

0.4252 

2. 372E-01 

2. 895E+00 

1 .85?9y 

41.6351 20.0000 1.719 

3.316E+00 

2. 241 E +00 1 . 479E+D0 8 

.44 3E-(I1 

2 . 800 

155.6 ! 

.401 3 

3. 6B4E-02 

9.460E-02 

0.3894 

2. 023E-01 

3.497E+00 

1.91480 

39.5046 20.0000 1.860 

3. 536E+00 

2.328E+00 1.519E+00 S 

.21 IE -01 

■! . noo 

142.6 1 

. 40 1 3 

2. 722E-02 

7.G24E-Q2 

0.3571 

1 .717E-01 

4. 229E+00 

1 .36463 

37.7787 20.0000 1.993 

3. 775E+00 

2.41 7E +00 1 . 562E+00 7 

.959E-U1 

3.200 

131.2 1 

.4012 

2. 024E-02 

6.169E-02 

0. 3280 

1 .452E-01 

5.1 12E+00 

2.00845 

36.3501 20.0000 2.119 

4.031E+00 

2.509E+00 1 . 607E+00 7 

.691E-01 

3.400 

120.7 1 

.4012 

1.514E-0? 

5.015E-02 

0.3018 

1 .226E-01 

6. 171E+00 

2.04707 

35.1484 20.0000 2.239 

4 . 904E + ni> 

2.601 E +00 1.G55E+00 7 

.41 0E-O1 

3.600 

111.3 1 

.4012 

1 . 140E-02 

4 . 056E-02 

0.2783 

1.035E-01 

7.431 £+00 

2.08122 

34.1245 20.0000 2.353 

4.593E+00 

2.693E+00 1 . 705E+00 7 

.120E-01 

3 . 600 

102.8 1 

.4012 

8.643E-03 

3. 362E-02 

0.2571 

8 . 743E-02 

8.924E+00 

2. M 149 

33.2431 20.0000 2.462' 

4 . 897E +00 

2. 785E+0U 1.750E+OO € 

. R23E-01 

4.000 

95.2 1 

.4012 

6. 599E-03 

2 . 773E-G2 

0.2380 

7. 397E-02 

1 . 068E+01 

2.13842 

32.4776 20.0000 2.566 

5.21 7E + 00 

2.876E+00 t . 814E+00 t 

. 523E-01 
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Computation of Thermally Perfect Properties of Oblique Shock Waves 


Table 6: Sample TPG tabular outputs for constant flow deflection angle 


(b) Strong shock wave solutions 


Thermally Perte:t Gas Properties Code I 
TPG, Version 5 . 1.0 I 

I 

NASA Langley Research rente' | 


Table of Thermally Perfe< » Compressible flow Properties 

Gaseous Mixture: Air: 4 -Species Mixture of N2. 02, Argon, and CO 2 

Database file name: [ Default database of Air Mixture ] 

Species Names are: 

HJ , 02 . Argon, C02 , 

Species Mass. Fractions are: 

0.75530 0.23145 0.01230 0.00040 

Species Mole Fractions are: 

0.78032 0.20946 t. 009 35 0.00026 

Mixture Properties: 

Molecular Weight - 28.3863 

Gas Cunstant = 2.87035E+02 J/(kq*Kj 

Polynomial Coefficients : cp/R =■ Suffi{A(i) + mi -2)} 

80.0 ■ T 1000.0 degrees h 

- 3. 7319821 0*00 5.4 1 %?rji[,-ni 3.46935860*00 3. 83537980-04 

-2.76531 300-06 8.18861420-09 -7.73686900-12 2.43487230-15 

1000. U T 6000.0 degrees h 

2.37638510*05 -1 . 2440527[i*C 3 5.12666900*00 -2.04008440-04 

6. 8321801 0- 08 - 1 .U553WJD- 1 1 6. 6450071 D- 16 0.00000000*00 

Total Temperature - 400.000 h 


ISENTROPIC FLOW PROPERTIES ---. 






OBLIQUE SHOCK WAVE 

HOW PRi. 

PERTIES 


2-dowrrst ream, 1-upstream) ----- 

M 

TO) 

Gamma 

P/Pt 

RHO/RHOt 

T/Tt 

q/Pt 

A/ A* 

V./3 + 

ShokAnql FloDfleC 

M2 

P2/P1 


RHCC/RHul 

T2/T1 Pt2/Pti 

0.000 

400.0 

1 .3949 

1.0006*00 

1 .OOOE+OO 

1 . 0000 

0. OOOE+OO 

inf ini te 

0.00000 

** Subsuni c : 

Nut Acini 

i cable 




0.200 

396.3 

1.3951 

9.7266-01 

9. 803E-01 

0.3922 

2. 714E-02 

2. 966E+00 

0.21784 

** Subsonic: 

Not Atiol i cable 





U. 400 

387.7 

1 . 3956 

8. 9596-01 

9.242E-0T 

0.9693 

1 . O0OE-O1 

1.531E+00 

0.43073 

** Subsonic: 

Not Aul 1 

icable 

... 



U. 600 

373.4 

1 . 3365 

7.844E-01 

6.403E-01 

0.9335 

1.372E-01 

1 . 189E+00 

0.63422 

** Subsonic: 

Not Applicable 





0.800 

355.0 

1 . 3975 

6.564E-01 

7. 397E-01 

0.8874 

2. 935E-01 

1 .033E+00 

0.82478 

+ * Subsonic: 

Not Aocli cable 





1 . 000 

333. 7 

1 . 3985 

5. 285E-01 

6.335E-01 

0.8343 

3.635E-01 

1 .0006+00 

1 .00000 

** Maximum At t ached -f 1 ow Shock 

Deflection Angle ■ 0 . 1300 

1.200 

310.3 

1 . 3994 

4.124E-01 

5. 305E-01 

0.7773 

4.155E-Q1 

1 .031E+00 

1.15866 

+ * Maximum Attached-flow Shock 

Deflection Angle « 3.953 

1 .400 

287.7 

1 .4001 

3. 141 E -01 

4.368E-01 

0.7192 

4.310E-01 

1 .115E+00 

1 . 30057 

+ * Maximum Attached- flow Shock 

De flection Angle ■ 9.444 

1 . 600 

264. 8 

1 .4006 

2.351E-01 

3. 551 E -01 

0.6620 

4.214E-01 

1 .251E+00 

1 .42631 

++ Maximum Attached-flow Shock 

Deflection Angle * 14.671 

1 . 800 

242.9 

1 .4010 

1 . 739E-01 

2. 863E-01 

0.6072 

3.946E-01 

1.440E+00 

1.53700 

** Maximum At t ached- f 1 ow Shock 

Deflection Analp » 19.201 

1 . 839 

238.7 

1.4010 

1 . 637E-01 

2. 743E-Q1 

0.5963 

3.873E-01 

1.483E+00 

1.55704 

64.9121 20.0000 

0.920 

3. 073E+00 

2 . 143E+00 

1.434E+00 8. 696E-C1 

2.000 

222. 3 

1.4012 

1 .2776-01 

2.297E-01 

0. 5553 

3.577E-01 

1 . 688E+00 

1 .63404 

74.3031 20.0000 

0. 728 

4.163E+ 

00 

2. 557E+00 

1 . 628E*QQ 7.555E-01 

2.200 

203. 3 

1 .4013 

9. 341 E -02 

1 .838E-01 

0.5083 

3. 168E-01 

2 . U05E+00 

1.71892 

77.5677 20.0000 

0.657 

5. 224E+00 

2. 882E+0D 

1 . 812E+00 6.517E-01 

2.400 

185.3 

1 .4013 

6.833E-02 

1.470E-01 

0.4648 

2.758E-01 

2.403E+00 

1.79311 

79.4141 20.0000 

0.613 

6.333E* 

00 

3.163E+00 

2. 002E+00 5.576E-01 

2.600 

1 70. 1 

1 .4013 

5.0086-02 

1.178E-01 

0.4252 

2.372E-01 

2. B95E+00 

1.85799 

00.6332 20.0000 

0.582 

7.519E* 

00 

3.41 IE *00 

2. 204E+00 4.735E-01 

2.800 

155.8 

1 .4013 

3.6846-02 

9.460E-02 

0. 3834 

2.023E-01 

3.497E+00 

1.91480 

81.5016 20.0000 

0 . 559 

8. 789E+00 

3. 633E*00 

2.41 9E *00 3.9S9E-01 

3. 000 

142.8 

1 .4013 

2.722E-0 2 

?. 624E-02 

U. 3571 

1.717E-01 

4 . 229E+00 

1.96463 

82. 1502 20.0000 

0.540 

1 .015E+01 

3. 832E+00 

2. 648E*00 3.367E-01 

3.200 

1 31 ,2 

1.4012 

2.024E-02 

6.169E-02 

0. 3280 

1.452E-01 

5.11 2E+00 

2. 00845 

82.6513 20.0000 

0.525 

1.159E+ 

01 

4.010E+00 

2.891E+00 2.830E-01 

3.400 

120. 7 

1 .4012 

1.514E-Q2 

5.015E-02 

0.3010 

1.226E-0I 

6.171E+00 

2.04707 

83.0484 20.0000 

0.519 

1.313E+ 

01 

4.1 70E+00 

3. 149E*00 2 . 377E-U1 

3. 600 

111.3 

1 .4012 

1 .1406-02 

4.096E-02 

0.2783 

1 . 035E-01 

7.431E+00 

2-08122 

83. 3695 20.0000 

o.5o: 

1.476E + 

01 

4.314E+00 

3.423E+00 1.998E-01 

3.800 

102.8 

1 .4012 

8. 643E-03 

3.362E-02 

0.25 7 ' 

8.743E-02 

8. 924E+00 

2.11149 

83.6336 20.0000 

0.493 

1.649E+ 

01 

4.443E+00 

3.711E+00 1.682E-01 

4.000 

95.2 

1 .4012 

6.593E-03 

2.7736-02 

0.2380 

397E-02 

1 .068E+01 

2.13842 

83.8536 2D.OUOO 

0 .486 

1.830E* 

01 

4.559E+00 

4.015E+00 1.419E-01 
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Computation of Thermally Perfect Properties of Oblique Shock Waves 


Table 7: Sample TPG output for a sequence of shock waves 


I Thermally Perfect Gas Proper tie? Code | 
I TPG , Version 3.1.0 I 
I I 
I NASA Langley Research center I 


Table nf Thermally Perfect Compressible Floe Properties 

Gaseous Mixture: Air: 4-Species Mixture of N2, 02. Argon, and C02 

database tile name: [ Default database of Air Mixture ] 

Species Names are: 

N2' , 02 , Argon, C 02 # 

Species Mass fractions are: 

0. 755.30 0.23140 I}. 01290 0.00040 

Species Me 1 e Frac Mon: are: 

0.78092 0.20946 0.00935 0.00026 


Mixture Properties: 

Molecular Weight = 23.9663 

Gas Constant 1 >.67035E*02 J/(kq*K) 
Polynomial Coefficients : cp/R = Sum} Ad >*t'‘u i -2)1 
Ki'JJj T ; 1000. Cl degrees K 

-3. 731 9821 [i*00 5.41967010-01 3.4693536D+00 

-2.76531 30b "06 3. 1 8861420 -09 -7.73636901' 12 

1000.0 T 6000.0 degrees h 
3 . 3769851 0 + 05 - 1 . 24405 270*03 5. 1 2G669I H«*(lO 

6.83218010-08 * 1 . 05535420-1 1 6.64500710-16 


83537980-04 

43407230-15 


2.041.106440-04 

0.00000000+00 


Total Temperature = 1500.000 F 


ISENTfiOPIC FLtW 

PROPERTIES 






OBLIQUE 

SHOCK WAVE 

FLOW PROPERTIES C 

.'■downst ream. 

1 "upstream) — > 

M 

T*>.) 

Gamma 

p/pt 

RHO/RHOt 

VTt 

g/Pt 

A/A'' 

V/a* 

Slink An g 

Fl o0f lee 

M2’ P2/P1 

RH02/RHOI 

T2/T1 Pt 2/P 1 1 

0.000 

1 500.0 

.3107 

i . on of + no 

1 . no OF +00 

i . nemo 

0. 010E+00 

inf ini te 

o.oouoo 


Subsonic: - 

Not Applicable — 



1 . 000 

1295.6 

. 313? 

5.422E-01 

6.277E-01 

0.8637 

3. 5756-01 

1 .000E+00 

1 . 00000 

** 

Maid mum Attached-flow Shock Reflection Angle = 0.000 

1.221 

1208.9 

. 2; 2 30 

4.077E-01 

5.059E-0’ 

Ij. 8059 

4. 035E-01 

1.042E+00 

1 . 19069 

71.1 36 

5. noun 

0.945 1 .408E+OO 

1 . 295E+00 1 

.fm&E + nn 9.956E-01 

7 . 11C0 

150.8 

.401 3 

2.0?3E-M4 

2. D66E-03 

0. 1005 

7. 1 35L-03 

1 . 233E+02 

2. 46 1 40 

1 1.827 

5. 000 Cl 

6.104 2. 235F+08 

1 . 749E + I20 1 

.27RE+00 9.491 E -01 



Total T e«ipereture = 1500 

000 K 









ISENTROPIf FLOW 

PROPERTIES 






OBLIQUE 

SHOCK WAVE 

FLOW 1 PROPERTIES ( 

-“downs tt earn. 

t -upstream) — ■ 

H 

Tt>:> 

Gamma 

P/Pt 

RHO/RHOt 

T/Tt 

q/Pt 

A/A*' 

V/a* 

ShokAngl 

F 1 L»[>f 1 e : 

M2 P2/P1 

RHLO/RHOI 

T2/T1 Pt >/Pt 1 

o.oon 

i son . n 

.3107 

1 . U00E +00 

1 . 000E +00 

1 . 0000 

0. OOOE+OO 

i nfirii te 

0.00000 

*«■ 

Subsonic: 

Not Applicable — 



0.945 

1315.1 

. 3179 

5 . 76SE-0I 

6. 579E-01 

0.8767 

3. 393E-01 

1 . 003E+00 

0.95158 


Subsonic: 

Not Appl i cab! e — 



1 . 000 

1295.6 

.3167 

5.422E Ml 

6.277E-01 

0.0637 

3.575E-01 

i.oorjE+oo 

1.00000 

** 

Maximum Attaclied-f 1 cvr Shock Deflection Angl 

e - fi.000 

1.231 

! 208 . 9 

. 3230 

4. 07 7 E -01 

5. 059E-01 

0.8059 

4. 085E-01 

1 . 042E+GCI 

1.19069 

71.1 361 

5.0000 

0.945 1 . 4Q3E+00 

1 . 295E+00 1 

. 038E+0C 9.956E-01 

6.104 

132. 7 

.4(ii 

4.895E-04 

3.81 IE -03 

0.1284 

1 . 278E-02 

6.789E+G1 

2.42635 

1.3.081 1 

5.0000 

5.398 2. 0 34 E + LIU 

1.643E+00 1 

. 238E +00 9.E.46E-01 
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Computation of Thermally Perfect Properties of Oblique Shock Waves 





Downstream Mach number M 


Computation of Thermally Perfect Properties of Oblique Shock Waves 



(a) Downstream Mach number M 2 and flow deflection angle 5 



(b) Density ratio p^/p; and temperature ratio T 2 /Tj 

Figure 3. Thermally and calorically perfect calculations of oblique shock wave properties within the 
calorically perfect temperature regime: air, 7^=400 K, a=30°. 
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Flow deflection angle 5 (deg) Temperature ratio T 2 / T, 



Static pressure ratio p 2 / p, F j ow Reflection angle 8 (deg) 






Computation of Thermally Perfect Properties of Oblique Shock Waves 



Figure 5. Effect of caloric imperfections on oblique shock wave properties within the calorically 
perfect temperature regime: air, T t - 400 K, a=30°. 



Figure 6. Effect of caloric imperfections on oblique shock wave properties within the calorically 
perfect temperature regime: air, T t = 400 K, a= 50°. 
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Downstream Mach number 


Computation of Thermally Perfect Properties of Oblique Shock Waves 



(a) Flow deflection angle 8 



(b) Downstream Mach number M 2 

Figure 7. Comparison of thermally and calorically perfect calculations for oblique shock waves in air 

at T,=1500 K for a=30° and 50°. 
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Temperature ratio T 2 / T, Pressure ratio p 2 / p, 







Total pressure ratio p u / Pl ,, Densjty ratjo Pj , p 


Computation of Thermally Perfect Properties of Oblique Shock Waves 


a = 50° 


v a = 30°- 


TPG 

0 equation 
CPG (a=30°) 
CPG (a=50°) 


Upstream Mach number M, 


(e) Density ratio P 2 /P 7 



~ TPG 

\ 

- 0 equation 

\ A 

CPG (o=30°) 

\ V 

CPG (o=50°) 


Upstream Mach number M 1 


(f) Total pressure ratio p t ^ P 
Figure 7. Concluded. 







V 


V 


V 


V 


r 


G = 30° 


— TPG 

- - 0 equation 

CPG (a=30°) 
CPG (o= 50°) 


J 1 1 1 U-J i_J i I JL I 




Temperature ratio T 2 / T, Pressure ratio p 2 / Pl 







Total pressure ratio 







(Mach,),, 

(Mach.) 


Upstream Mach number M 1 


(a) Flow deflection angle 6 



123456789 10 


Upstream Mach number M, 

(b) Downstream Mach number Afg 

Figure 9. Effect of caloric imperfections on oblique shock wave properties: comparisons between TPG 
and NACA Rep. 1135 for two total temperatures and two shock wave angles 






Computation of Thermally Perfect Properties 

of Oblique Shock Waves 

TPG 


© eqn. 

/y 

// 

// 


// 

// 

// 

\ 

h 

it 

ft 

\ T t = 3000 K 

it 

it 

t a = 50°^^^^^^ 


j/ /^a^O 0 

\ 


> T,= 1500 K 


Upstream Mach number M 1 

(e) Density ratio p^/p i 


a = 30° 



T = 1500 K 



Upstream Mach number M t 


T. = 3000 K 


(f) Total pressure ratio Pt^PtJ 
Figure 9. Concluded . 






Lateral distance from Wedge Leading Edge 


Computation of Thermally Perfect Properties of Oblique Shock Waves 



TPG shock wave 


Wedge angle, 5 = 10° 


Longitudinal distance from Wedge Leading Edge 

Figure 10. GASP 2-D Euler solution pressure ratio contours and TPG oblique shock wave for 
thermally perfect air: Afj=4.8, 7\=1670 K over a 10° wedge. 


GASP 2j*D Euler at z=2.6 
TPG Solution 


Longitudinal distance from Wedge Leading Edge 


(a) Mach number M 

Figure 11. GASP 2-D Euler and TPG oblique shock wave calculations for thermally perfect air: 
Mi= 4.8, T;=1670 K over a 10° wedge: at z=2.6 units above the wedge leading edge. 











Lateral distance from Wedge Leading Edge 







Pressure ratio p / p, Mach number 


Computation of Thermally Perfect Properties of Oblique Shock Waves 



Longitudinal distance from Wedge Leading Edge 

(a) Mach number M 

5 
4 
3 
2 


1 

5 6 7 8 9 

Longitudinal distance from Wedge Leading Edge 

(b) Static pressure ratio p/p 2 

Figure 13. GASP 2-D Euler and TPG oblique shock wave calculations for 20% Steam/80% C0 2} by 
mass: M 2 =3.5, T 2 = 1200 K over a 25° wedge: at z=5.0 units above the wedge leading edge. 

43 









REPORT DOCUMENTATION PAGE 


Form Approved 
OMBNo. 0704-0188 


Public reporting burden for this collection of information s estimated to averaop i hour nor r«rwi«o i-.*, ■ 

gathering and mamtanmg the data needed and comofet^a and review™ ^ * reviewing instructions, searching existing data sources. 

1. AGENCY USE ONLY {Leave blank) I 2. REPORT DATE ^“REPORTTYPE AND DATPS rnucocn 


August 1996 


4. TTTLE AND SUBTITLE 


3. REPORT TYPE AND DATES COVERED 
Contractor Report 


Computation of Thermally Perfect Properties of Oblique Shock Waves 


5. FUNDING NUMBERS 


C NASI -19000 


6. AUTHOR(S) 

Kenneth E. Tatum 


WU 505-70-62-01 


7. PERFORMING ORGANIZATION NAM£{S) AND ADDRESSEES) 

Lockheed Martin Engineering & Sciences Company 
Hampton, VA 23666 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESSEES) 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23681-0001 


10. SPONSORING /MONITORING 
AGENCY REPORT NUMBER 

NASA CR-4749 


11. SUPPLEMEMTARY NOTES 


Langley Technical Monitor R. Wayne Guy 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 

I Unclassified - Unlimited 


12b. DISTRIBUTION CODE 


Subject Category 02 


13. ABSTRACT (Maximum 200 words) — — 

A set of compressible flow relations describing flow properties across oblique shock waves, derived for a 
thermally perfect, calorically imperfect gas, is applied within the existing thermally perfect gas (TPG) computer 
code. The relations are based upon a value of cp expressed as a polynomial function of temperature The 
updated code produces tables of compressible flow properties of oblique shock waves, as well as the original 
properties of normal shock waves and basic isentropic flow, in a format similar to the tables for normal shock 
waves found in NACA Rep. 1135. The code results are validated in both the calorically perfect and the 

imperiect ’ thermal| y P erfect temperature regimes through comparisons with the theoretical methods of 
NACA Rep. 1135, and with a state-of-the-art computational fluid dynamics code. The advantages of the TPG 
code for oblique shock wave calculations, as well as for the properties of isentropic flow and normal shock 

waves, are its ease of use, and its applicability to any type of gas (monatomic, diatomic, triatomic, polyatomic or 
any specified mixture thereof). 


14. SUBJECT TERMS 

Thermally perfect gas, calorically imperfect gas, thermally perfect gas computer 
code, compressible flow relations, oblique shock waves, imperfect gas effects. 


115. NUMBER OF PAGES 


j 16. PRICE CODE 


1 17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


120. LIMITATION OF ABSTRACT 


NSN 7540-01 -280-5500 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. 239-18 
298-102 



